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The Einstein- Vlasov-Maxwell(EVM) system can be viewed as a rela- 
tivists generalization of the Vlasov-Poisson(VP) system. As it is proved 
below, one of nice property obeys by the first system is that the strong 
energy condition holds and this allows to conclude that the above system 
is physically viable. We show in this paper that in the context of spherical 
symmetry, solutions of the perturbed (EVM) system by 7 := 1/c 2 , c being 
the speed of light, exist and converge uniformly in the L°°-norm, as c goes 
to infinity on compact time intervals to solutions of the non-relativistic 
(VP) system. 

1 Introduction 

The classical (VP) system models the time evolution of collisionless particles 
in the Newtonian dynamic setting, particles could be for instance atoms and 
molecules in neutral gas or electrons and ions in a plasma. In stellar dynamics, 
particles are either stars in galaxy or galaxies in a cluster of galaxies pQ. The 
global behaviour of solutions of the above system now is well understood (see 



EH- [TTj . |T2| . [IT)] . 0). In \T^, the authors use an estimate to prove that 



the Newtonian limit of the spherically symmetric Vlasov-Einstein system is the 
classical (VP) system. This result is extended by Rendall in ^1] to the general 
asymptotically flat Einstein- Vlasov system. The above investigation concerns 
the case where the charge of particles is small to be neglected. 

We are inspired by what is done in and we want to extend this result 
considering in this paper the (EVM) system that models the time evolution of 
self-gravitating collisionless charged particles in the general relativity setting. 
Firstly we discuss and obtain that the above system satisfies the strong energy 
condition, and then the system is physically viable. Secondly, we perturb the 
(EVM) system by a parameter 7 = 1/c 2 and together with the assumption of 
spherical symmetry, we show using an estimate that for 7 small, the Cauchy 
problem associated with the obtained system admits a unique regular solution, 



symmetric Einstein- Vlasov-Maxwell 
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with time existence interval independent of 7. Once this result is obtained we 
deduce that if c goes to infinity then the above solution converges uniformly to a 
solution of the classical (VP) system. To do so, we consider the new constraint 
equations and discuss the existence of solutions satisfying the constraints as 
we did in [7|- This gives rise to a new mathematical feature since the present 
equations contain two parameters 7 and q (q being the charge of particle), 
rather than one as it is the case where 7=1. We obtain that for given initial 
datum to the distribution function, solutions of the constraints exist and depend 
smoothly(i.e C°°) on parameter £ := (7,9), for £ small, and this allows us to 
construct the set of initial data so that the initial datum for the metric function 
A 7 is bounded in the L°°-norm. So, with these initial data, we undertake the 
Cauchy problem for the modified (EVM) system and establish as we did in [S], 
the local existence and uniqueness theorem and continuation criterion. The 
interest of this paper lies on the fact that the estimates we use are complicated 
to establish. In fact, contrary to the uncharged case, due to the presence of 
metric function A 7 in the both side of equations, we have to estimate both the 
supremum of momenta on the support of the distribution function / 7 and the 
supremum of e 2X ~< . We are not aware that the above is already done. 

It is appropriate at this point to put our investigation in the context of 
general relativity. The classical limit of the relativistic Vlasov-Maxwell system 
is studied by Jack Schaeffer [T7j. Using and estimate, the author proves that 
this limit is the classical (VP) system. Also, the same result is established by 
Asano for the classical Vlasov-Maxwell system. The recent result in this field 
is that of S. Calogero and H.Lee 0- They consider the relativistic Nordstrom- 
Vlasov system and establish that its non-relativistic limit is the classical (VP) 
system. 

The paper is organized as follows. In Sect. 2, we recall the classical (VP) 
system. In Sect. 3, we introduce the general formulation of the (EVM) system 
and we discuss the energy conditions. In Sect. 4, we introduce the perturbed 
spherically symmetric (EVM) system and state the main results of this paper. 

2 The classical(non-relativistic) (VP) system 

In the Newtonian dynamic, the time evolution of a set of collisionless particles 
is governed by the following equations known as the (VP) system: 



Here (|2.1|) is the Vlasov equation for the unknown /, / being the distribution 
function that measures the probability density to find a particle (star) at time 
t with position x and with momentum p, where < > 0, £ I 2 . Note that 

/ is defined on the mass shell. (|2.2J) is the Poisson equation for the unknown 



d t f + v-V £ f-VzU-Vpf = 



(2.1) 




(2.2) 
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U = U(t,x) that measures the Newtonian potential generated by stars. If 
r\ = — 1 then (|2.1() and l|2.2l) model the plasma physics case. In what follows, we 
consider the case where r\ = 1. 

3 The (EVM) system 

As we said before, we take fast moving collisionless particles with unit mass 
and charge q. The gravitational constant and the speed of light are taken to 
be equal to unity. The basic spacetime is (R 4 ,g), with g a Lorentzian metric 
with signature (—,+,+,+). In the sequel, we assume that Greek indices run 
from to 3 and Latin indices from 1 to 3. We adopt the Einstein summation 
convention. The (EVM) system reads in local coordinates at a given point 
(x<*) = (t,x): 

Ra0 - \g a( iR = MT a0 (f) + r a p(F)) (3.1) 
df p 1 df { a g 4 df 

m + V^" {TpJp + qp Fp ] W = (3 ' 2) 

V a F af3 = J p {f)- V a F M + V B F^ a + V 7 F aB = 0, (3.3) 

with: 

T a p{f) = - I PaPpUJp] T a p 

J (f)(x)=q I p^f(x,p)LJ p ; LO t 

JR 3 

where T"^ denote the Christofcll symbols. In the above, / stands for the distri- 
bution function of the charged particles defined on the mass shell: 

g a p P a P P = -1, 

F stands for the electromagnetic field created by the charged particles. Here 
l|3.1|l is the Einstein equations for the metric tensor g = (g a g) with sources 
generated by both / and F, that appear in the stress-energy tensor T a g + r a g. 
Equation (|3.2|) is the Vlasov equation for / and (|3.3|l are the Maxwell equations. 

3.1 The energy conditions 

Physically, the quantity (T a g + Ta^V^V 13 represents the energy density of 
charged particles obtained by an observer whose 4- velocity is (V a ). So, for 
any physically viable theory, this quantity is nonnegative for every timelike vec- 
tor (V a ), and the above assumption is known as the weak energy condition. 
We are going to show that in fact, the dominant energy condition holds, i.e: 
{T a 3+T a0 )V a W 13 > 0, for all future-pointing timelike vectors (V) a and (W a ). 
This implies the weak energy condition. Note also that the above definition 



H/2 dp 1 dp 2 dp 3 

9 I ' , Po = 9ooP , 

Po 
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of dominant energy condition is equivalent to that given by Hawking and El- 
lis in [ 5 , p. 91]. Once this is clarified, we show that since the Maxwell tensor 
is traceless, and the strong energy condition that is R a pV a VP > 0, for every 
timelike vector (V a ), holds for the Einstein- Vlasov system, the same is true in 
our context. 

Lemma 3.1 The following assertions are equivalent: 

1) For any two future-pointing timelike vectors (V a ), (W a ), one has: 
(T af3 + T aP )V a WP > 

2) For every timelike vector (V a ), one has: (T a p + T a p)V a V^ > 0, and 
{T a p + T a p)V^ is a non- spacelike vector. 

Proof: We first prove that 2) implies 1). If (V a ) and (W a ) are two future- 
pointing timelike vectors then 2) implies that (T a p + T a p)V^ is non-spacelike. 
This means by definition that it is either timelike or null. Using once again 2), 
its contraction with (V") is positive. Thus, it is in fact past-pointing timelike or 
null, its opposite is future-pointing timelike and then its contraction with (W a ) 
is positive. This comes from the fact that a contraction of two future-pointing 
timelike vectors is negative or null. So 1) is proved. 

Now conversely, suppose that 1) holds and let us prove 2). Let (V a ) be 
a timelike vector. If V° > 0, then (V a ) is future-pointing timelike and the 
first condition in 2) holds by taking V = W. If V° < 0, then (-V°,-V l ) is 
future-pointing timelike and we can conclude as we made before for the first 
condition in 2). For the second part, suppose that (V a ) is future-pointing 
timelike and define P a := {T a p + r a p)V^. Condition 1) implies that (P a ) 
satisfies P a W a > for every future-pointing timelike vector (W a ). We aim 
to show that (P a ) is non-spacelike. To do this, let us assume that (P a ) is 
spacelike, and get a contradiction. Set L := P a P a . By assumption L > 0. Let 
(T a ) be a future-pointing timelike vector orthogonal to (P a ) with T a T a = —L 
and T° > P°(for the construction of vector (T a ), one can take for instance in 

normal coordinates: T° = Jf^P 4 ) 2 , T l = P°P l /T°). Set W a = 2T a - P a . 

Then W a W a = -3L < 0, W° = 2T° - P° > and (W a ) is a future-pointing 
timelike vector, and W a P a = —L < 0. This is the desired contradiction. Now 
if (V a ) is past-pointing timelike, then (— V a ) is future-pointing and follow the 
first step of the proof in which P a is replaced by — P a = {T a fj + T a p){— V 13 ). 
Analogously, we are led to (—P a ) is non-spacelike and so is (P a ). In conclusion, 
the second part of condition 2) holds, for every timelike vector (V a ) and the 
proof is complete. 

Proposition 3.1 1) For every two future-pointing vectors (V a ), (W a ), one 
has: 

T a pV a W + T aP V a W p > (3.4) 
2) For every timelike vector (V a ), one has: 

R aP V a V p > 0. 
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Proof: Consider part 1) of the above proposition. Since Penrose and Rindlcr 
state the dominant energy condition for the Maxwell equations in writing the 
Maxwell tensor r a p as a quadratic form of spinor fields in JU| , the second term 
in the left hand side of l|3.4|) is nonnegative and we just need to establish the 
same result for the first term in the left hand side of (|3.4|) . Let (V a ), (W a ) 
be two future-pointing timelike vectors. Taking the first term in the left hand 
side of H3.4|) . we obtain, since / > 0, (p a ) is future-pointing timelike vector and 
— Po > 0, one has: 

TaflVW' = f ( Pa V°)(ppWP)f | g \i dpldp2dp3 > 
Jm.3 -P0 

So, (|3.4[1 holds as announced. Now concerning part 2) of the above Proposition, 
the contraction of (|3.1[1 gives, since g a/3 T a p = 0: 

R = -8vrT 

where T := g a/3 T a p. Insertion of the above in (|3.1|) yields: 

Rap = -iirTg a p + &n(T a p + Tap). 
Next, let (V a ) be a timelike vector. Then 

R a pV a V f3 = ^T af) V a V p + 8tt (r aP - \Tg a ^j V a V? (3.5) 

Since the Maxwell tensor satisfies the dominant energy condition and then the 
weak energy condition, we can deduce that the first term in the right hand 
side of (|3.5J) is nonnegative. Also, the strong energy condition holds for the 
Einstein- Vlasov system (for more details one can refer to [|15|. p.37-38]). The 
latter shows that the second term in the right hand side of (|3.5() is nonnegative 
and the strong energy condition holds in our context. 



4 The perturbed spherically symmetric (EVM) 
system and the main results 

We consider the Lorentzian spacetime (R , g 7 ) where g 7 is obtained by scaling g 
with 7 = 1/c 2 in the spatial part. With the assumption of spherical symmetry, 
we take g 7 of the following form: 

ds 2 = -e^dt 2 + je 2X dr 2 + r 2 j(d9 2 + sin 2 0dip 2 ) 

where /i = n(t,r); X = X(t,r); t eR;r e [0, +oo[; 6 e [0,tt]; ip £ [0,2tt]. 

Next, the assumption of spherical symmetry allows us to deduce the following 
perturbed (EVM) system in the (t, x, v) coordinates (see |H]); where 

v l := p* + (e A - 1)^- : 
r r 
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e~ 2A (2rA' - 1) + 1 = 87r 7 r 2 p 
e" 2A (2r j u / - 1) + 1 = 8TTj 2 r 2 p 



(4.1) 
(4.2) 



ot J\ + -yii 2 \ 7 r J r ov 

(4.3) 

— (r 2 e A e) = gr 2 7 3/2 e A M (4.4) 
where A' = f ; A = § and: 

p(t,x):= [ f(t, x, v)y/l + 7 v 2 dv + le 2Mt ^e 2 {t, x) (4.5) 

JR3 2 

p(t,x):=[ (^\\ { t^v)^=-\e 2 ^e 2 {t,i) (4.6) 

M(t,x) := { f(t,x,v)dv. (4.7) 



Next, it is well known that in the context of spherical symmetry, the (VP) 
system reduces to: 

d t f + vdzf-K(t,x)-d v f = (4.8) 



K(t,x) = \- f M(t,y)dy. (4.9) 
r r J\y\<r 

In the sequel we denote by (EVM 7 ) the system (|4.1|l . (|4.2|) . (|4.3|l and l|4.4|) . 

o o o 

It is prescribed on (£yM 7 ) the initial data A(0) = A and /(0) = /, where / 
is nonnegative, spherically symmetric function i.e invariant when one applies 

o 

any rotation of R 3 on both the variables i and v, A being a given function. 

o 

Concerning the system (14.811 and l|4.9|l . only the initial condition /(0) = / is 
needed. Also, we are interesting with the asymptotic flatness of spacetime and 
a regular center at r = 0. So, we have the following boundary conditions: 

lim X(t,r)= lim fi(t, r) = A(0, r) = 0, t > 0. (4.10) 

r — »+oo r — >-foo 

Now, the assumption that the electric field E defined by E(t,x) — e(t,r)-, is 
spherically symmetric and the fact that at the spatial infinity there is no charged 
particle lead to the following boundary condition (for more details, see 0): 

lim e(i, r) = e(0, r) = 0, t > 0. (4.11) 
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Before stating the main result of this section, we recall the following result whose 
proof is made by straightforward calculation (see |§]). The concept of regularity 
of solutions we use is the same as in [5] . 

Lemma 4.1 Take 7 = 1. Let (/, A, fi, e) be a regular solution of EVM\ satis- 
fying \4.lU\j and \4-H\l on some interval I . Then for every a£l \ {0}, 

f f a (t, x, v) := a 2 f(at, ax, v); A a (t, r) := \(at, ar) 
\na(t, r) := n(at, ar); e a (t, r) := ae(at, ar) 

defines another regular solution of this system on the interval a -1 /. 

Another thing to discuss is the existence of initial data satisfying the constraints. 
The constraint equations are obtained by setting t = in (|4.1|l . i|4.2[l and l|4.4|) . 
As we said in [7J, these constraints reduce to 

e _2A (2rA - 1) + 1 = 87r 7 r 2 p (4.12) 

4(r 2 e S e)= 9 r 2 7 3/ VM (4.13) 
dr 

where 

e(r) := e(0, r); p(r) := p(0, r); M(r) := M(0, r) 
If 7 G] 0,1] and 

8tt / s 2 ds / °f(s,v)y/i + v 2 dv < r, r > 0, (4.14) 

then 



i) 



87T7 / s 2 ds / /(s, u) yl + r yv 2 dv < r, r > 0. 



So, under the assumption l|4.14|l . we can apply the results obtained in [7] to 
H4.12jl and (|4.13j) , and deduce the following 

o 

Proposition 4.1 Let f G C°°(R 6 ) be nonnegative, compactly supported, spher- 
ically symmetric and satisfies \4-.14\j . Take 7 G]0, 1]. Then, for Q = (7,9) 
small enough, the equations \4-TSfy and have a unique global solution 

(A 7 ,e 7 ) S (C°°([0, +oo[)) 2 such that A 7 (0) = e 7 (0) = 0. Moreover, this solution 
depends smoothly (i.e C°° ) on the parameter £. 

Next, with the help of Proposition 4.1, for 7 G]0, 1] and A > 0, one can define 
the set of initial data as: 

D := {(/,A 7 ,e 7 ) e C°°(IR 6 ) x (C^^+ooQ) 2 , / > 0, spherically symmetric, 
and satisfies (14.14(1 and (A 7 ,e 7 ) 

o 

is a regular solution of the constraints with |j A 7 ||l^< A}. 
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We now study the Cauchy problem given by (EVM 7 ), the initial data, 14.1L)f> 
and (EU). 

Theorem 4.1 There exists T > and a continuous function u : [0, T[— > M+ 

° ° o 

sucft that for 7 s]0, 1] and (/, A 7 , e 7 ) € D, i/ie system (EVM^) has a unique reg- 
ular solution (/ 7 , A 7 , /x 7 , e 7 ) on £/ie interval [0, T[ w«<A initial data (/, A 7 , /z 7 , e 7 ) 
and 

/ 7 (*,x,u) = 0, \v\>u(t), ieR 3 ,ie[0,T[, 
where /x 7 := /^ 7 (0, r). 

00 

Proof: Take 7 G]0, 1] and (/,A 7 ,e 7 ) G D. For 7 = 1, it is shown in [S] 
that the system (£VMi) admits a unique regular local solution (/, A, fi, e) de- 
fined on a maximal existence interval [0, Ti[, with initial data (f,Xi,fi 1 ,e\), 
and using Lemma 4.1, one deduces that (/ 7 := 7 3 / 2 /(7 1 / 2 ., 7 1 / 2 ., .), A 7 := 

A(7 1 / 2 .,7 1 / 2 .),Ai 7 := M7 1/2 -,7 1/2 -),e 7 := 7 1/2 e( 7 1 / 2 ., 7 1/2 .)) solves (£FM 7 ) 
on [0,T 7 [, where T 7 := 7~ 1/2 Ti. Set for every t G [0,T 7 [ 

U y (t) := sup{| « I |(x,u) G supp/ 7 (t)} 
Q 7 (f) :=sup{e 2A - (t ' r) |r > 0}. 

We now look for an estimate for U 7 and <5 7 . Let s t— > (X~(s, t, i, u), V^s, i, x, v)) 
be the solution of characteristic system 

V = - I A 7 — + e^~ A ^-u' Jl + 71; 2 - qe^ +A ^e, 
\ r 7 ' 

with Xy(<;, f, x, v) = x and Vj(t,t,x,v) = v. It is well known that the solution 

o 

of 14. 3|) with initial datum / is given by 

o 

f y (t,x,v) = f(Xj(0,t,x,v),Vy(0,t,x,v)) 

o 

and || f-y(t) ||j,oo=|| / ||l=c, t > 0. Along a characteristic, we obtain since 
^ 7 — A 7 < 0: I x \<\ v I / \J"{v 2 < -j= and one deduces from the above by 
integration on [0,t] that r :=| i |< r + where 

o 

r := sup{| x I I (x, v) G supp/} 
The integration of (|4.4|l yields: 

e 7 (t,r) = A 7 3/2 e "A T (t,r) /" s V^ s) M 7 (M)ds (4.15) 
r 2 7 
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where M 7 is deduced from (|4.7II . replacing / by / 7 . Now, since — A 7 < 0, 
distinguishing the cases r < ro + A= and r > tq + -4= , one obtains the following 

estimate for e 7 , since 7 g]0, 1]: 

|e 7 (i,r) \<C{l + r + t)Q^ 2 {t)U^{t), r>0 (4.16) 

o 

wl : C = C(q, II / ||i°o) is a constant. Using (|4.16|l . we deduce 

(e 2A -e 2 )(i, r) < C(l + r Q + t) 2 Q 3 J 2 (t)U®(t), r > (4.17) 

We use H4.17|l and estimates obtained in the proof of Theorem 1 in JSj to deduce 
and estimates for p 7 and p 1 : 



II P-y(t) Hi- < CU*{t)^\ + U 2 {t) + C(l + r + i) 2 g 7 /2 (t)^(i), 
lb 7 W <Cmin|-^C/ 7 4 (t),L/ 7 5 (t)|+C(l + r + t) 2 Q 7 / 2 WC/ 7 6 (t). 

Here p 7 and p 7 are deduced from 14.5|l and 1)4. 6[) replacing /, A and e by / 7 , A 7 
and e 7 respectively. We now introduce the mass function given by: 



m 7 (t, r) :— An s p 7 (t, s)ds = / p~ l {t 1 y)dy 

Jo J\v\<r 

and with the help of the estimate of p 7 , one obtains, after distinguishing the 
cases r < r + -7= and r > r + -4=, and using the relation 

v t v T 

(e A ^e 7 )(t,r) = ( 7 ° + *^ ) (e A -e 7 )(t, r + t/^y), r G [r + t/V7>+°°[ 
Now, one has: 

W7 ^' ^ < 4tt / s Pl {t,s)ds 
r Jo 

<c(r a + -^j U*{t)^l + f/ 7 2 (t) + C (r + -^=) (1 + r + i) 2 Q 7 /2 (7 7 6 W 

< ^ (r + -1) f/ 7 3 (t) V / TT^) + C (r + J=) (1 + r + i) 2 Q 7 / 2 ^(i) 

The integration of (|4.1|l on [0,r] yields e _2A T as a function of m 7 , and the 
insertion of this relation in l|4.2|) gives ji' . So multiplying the obtained equation 
by e^~ x f yields: 

(e^- A >;)(i,r) = e^ +A - ~ + 47r 7 2 rp 7 (t, r)) . 
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Since /z 7 + A 7 < 1, one obtains the following estimate, after distinguishing the 
cases r < tq + tj and r > r$ + tj ^/j for the second term in the right hand 
side of the expression above: 



7 7 



< £7(1 + r Q + t)U*(t)yJl + U*(t) + £7(1 + r a + t) 3 C£/ 2 (t)U% (t) 
+ £7(1 + r + i)C/ 7 5 (;) + (1 + r + i) 3 Q 7 (t)t/ 7 6 (i). 

(4.18) 

Another part of the Einstein equations is (see [5]): 

Xy(t,r) = -4nyre' 1 *' +> *ky(t,r), 

where 

f x .v 

k 7 (t, r) — k 7 (t, x) := / f 7 (t,x,v)dv, 

from which one deduces the following estimate for A 7 : 

I A 7 (i,r) |< £7(1 + r + t)U*(t) (4.19) 
Inserting l|4.16fl . I|4.18|) and (|4.19() in the characteristic system, one has: 
I v | < £7(1 + r Q + t)U*{t) + £7(1 + r Q + t)U 3 (t)(l + U*{t)) 
+ £7(1 + r + t) 3 QV 2 (t)U*(t)y/l + U*(t) 



+ £7(1 + r + t)U*(t)yJl + U*{t) 

+ £7(1 + r + t) 3 Q*(t)U%(t)yJl + U*(t) 
+ C(l + r +t)Q 1 / 2 (t)U 3 (t). 



Since 



Q 7 /2 W < yJl + Q-y® < 1 + Q 7 (t), 

q 7 /2 (*) < (i + g 7 w) 3 < (i + q 7 (*) + ^(*)) 3 

one deduces from the above that 

| 6 |< £7(1 + r + i) 3 (l + Q 7 (i) + C/ 7 (i)) n - 
So, the integration of the above inequality on [0, t] gives: 



U 7 (t) <U + C (1 + r + s) 3 (l + Q 7 (s) + C/ 7 (s)) n ds, (4.20) 
Jo 

where 

o 

U := sup{| « | | (5, u) G supp/}. 
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We recall that in the above, the constant C depends on q and || / H^oc and does 
not depend neither on t and nor on 7. (|4.2U|) shows that we need an estimate 
for Q 7 (t). We will proceed exactly like we did in [5] when showing that the 
sequence of iterates is bounded in the L°°-norm. So, we differentiate e~ A ^ w.r.t 
t and obtain: 



d 



,2A T (t,r) 



< 2jQ 2 Jt) 



m 7 (t, r) 



We see from the above inequality that an estimate for is needed. Using the 
Gauss theorem after insertion of the Vlasov equation, we deduce: 



TYl™ .. \ . 1 



\y\<r JR 3 I V 



\y\<r JR 3 



\y\<r 



dvdy 



-a-,./ V- v 



1 \y\ 



(4.21) 



Q 
r 

2n 

r 



\y\<r JR 3 



y.v 



--f-ydvdy 



- / , s 2 | ( e 2 



Another part of the Maxwell equations is given by: 

l( e ^e 7 ) = -V/V^ 7 



where 



AT 7 (t,r) = N 7 (t,x) 



Ir 3 \/r+ ■yv 2 

So, the last term in the right hand side of l|4.21(l yields 



/ 7 (i, x, v)dv. 



2tt 
r 



dt 



(e 2A ^e 2 )ds = -47T 7 3 / 2 ^ / se^+^e^ds, 



and since 



JV 7 (i,r) |< 



fC(r /V7 + t/7)^ 7 3 W 
\C(r +i/V7)C/ 7 4 W, 



we deduce the following estimate: 



2tt 
r 



S 2 |(e 2 ^e 2 )d S 



<C(l + r + i) 3 Qy 2 (t)^ 6 W. 



(4.21a) 



Next, we give an estimate for the first term in the right hand side of (|4.21|) . 
after distinguishing the cases r < rp + t/y/y and r > tq + t/^/y: 



>k y \<C(l+r +t)U*(t) 



(4.21b) 
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Next, from 1)4. lfl . we deduce the following 

e^-M; = e^ +x -< (-7^ + 4vr 7 rp 7 ) 

and then 

e^- x -<{^ - A' 7 ) = e^ +A ^ (27^ + 4vr7 2 rp 7 - 47r 7 rp 7 ) , 
and using the above , one has: 

I e^" A - (/4 - A;) I (t, r) < C(l + r + t)U^l + U*(t) 

+ C{\ + r +tfQ^ 2 {t)U^{t) 
+ C(l + r +t)U*{t). 

So, with the help of this, and denoting by A the second term in the right hand 
side of (|4.21() . one deduces the following estimate for A: 

t ^ 2 



\A\<C(l + r + t)[r + —) U?{t) 

+ <7(l + r + t)(Vo + -j=) U%(t)yjl + U*(t) 

+ C(l + r +t) 3 (V + -j=) C£/ 2 U?(t). (4.21c) 

Using H4.19|l . we obtain the following estimate for the third term, say B, in the 
right hand side of (|4.21(l : 

\B\<C(l + r +t)(ro + -^j U*(t)y/l + U*(t) 

+ C(l + r +t)(r + ^j C7 7 °(t). (4.21d) 

Using H4.18|l . we deduce an estimate for the fourth term, say C, in the right 
hand side of ijOT)) : 

\C\< C 7 (l + r + t) (r + Ul(t)yjl + U*(t) 

+ C 7 (l + r + tf (r + Qy 2 {t)U}?(t) 
+ C 7 (l + r + t) (ro + ^J C/ 7 9 (i) 

+ C 7 (l + r + i) 3 (»"0 + -^) Q'(*)^°(*)- ( 4 - 21e ) 
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Using H4.15fl . the fifth term, say D in the right hand side of (|4.21|) can be 

estimated by: 

\D\<C 1 3 ' 2 (r + ±^ Q^iW^t). (4.21f) 

So, taking into account (4.21a), (4.21b), (4.21c), (4.21d), (4.21e), (4.21f), one 
obtains the following estimate: 



d 



,2A 7 (t,r) 



<C(l + r + t) 3 Q^ 2 (t)U^t) 
+ C(l + r Q +t)Q 2 (t)U 3 (t) 



+ C{1 + r Q + t) 3 Q 2 (t)U»(t)^Jl + U 2 (t) 

+ C(l+r +t) 5 Q^ 2 (t)U} / 2 (t) 
+ C(l + r +t) 3 Q 2 (t)U? ( 1 (t) 

C(l +r + t) 3 Q 2 (t)U*(t)^/l + U 2 (t) 
+ C(l + r +t) 3 Q 2 (t)U™(t) 



+ C(l + r + t) 3 Q 2 (t)U^t)^l + U 2 (t) 

+ C{l + r +tfQy 2 {t)U™{t) 
+ C(l + r +t) 3 Q 2 (t)U»(t) 
+ C(l + r +t) 5 Q*(t)U™(t) 
+ C(l + r Q +t) i Q^ 2 (t)U^(t), 



thus, 



d 

8? 



,2A T (t,r) 



\17 



< C(l + r + t) 5 (l + Q 7 (t) + U 7 (t)Y 
So, the integration of the above inequality on [0,t] yields: 

g 2A 7 (f,r) < e 2A 7 (r) + Q f (1 + f(j + fi)5(1 + g^ fl) + f/^))^ 



and bearing in mind that (/, A 7 , e 7 ) S I?, one has an estimate for Q 7 : 



(t)<C + C (1 + r + s) 5 (l + Q 7 (s) + C/ 7 (s)) 17 ds, 



(4.22) 



where C is a constant which depends on /, q and A, and not on 7 and t. Next, 
adding l|4.20|l and (|4.22|) . one has: 

Q 7 (t) + l/ 7 (t) <U + C + C [ (1 + r + s) 5 (l + Q 7 (s) + C/ 7 (s)) 17 ds, 
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for t £ [0,T 7 [. Let u : [0,T[ 
integral equation: 



be the maximal solution of the following 



u(t) = U a + C + C (1 + r + s) b (l + u(s)) u ds. 



Then 



U y (t) < Q^(i) + fXy(i) <u(t), t£ [0,T 7 [n[0,T[ 

and with the continuation criterion proved in |S] , one concludes that U 1 (t) < u(t) 
for t £ [0, T[ and this ends the proof of Theorem 4.1. 

We now state and prove the essential result of this paper. This is concerned 
with the convergence of a solution (/ 7 , A 7 , fi~, e 7 ) of (EVM^) to a solution / 
for the (VP) system given by l|4.8|l and Q4.9p . as 7 tends to 0. 

Theorem 4.2 Let < T < 00 6e smc/i ttai /or every 7 G [0, 1[ and 

o 

(/,A 7 ,e 7 ) € D, f/ie solution (/ 7 , A 7 , /it 7 , e 7 ) of (EVM 7 ) exists on [0,T[ and 

f 7 (t,x,v) = 0, I t> |> u(i), ieR 3 ,te [0,T[, 
where u : [0,T[— > K + is a continuous function. Let f £ C 1 ([0, +oo[xR 6 ) 6e 

o 

£/ie solution of (VP) with /(0) = /. TTien /or every T' G]0,T[ i/iere exists a 
constant C > swc/i £/ia£ /or any 7 G]0, 1] i/ie following estimate holds: 



II IU» 
+ II e 7 W IU« 

/or every t £ [0, T'] . 



A 7 (t) 



~d 2 /i 7 (t) - JT(t) 

7 



<C 7 , 



A 7 (i) 



(4.23) 



Proof: Take 0<T' <T. Let Ci be the upper bound of u on [0, T']. Then one 
obtains, using the estimates above: 

U,(t), Q 7 (i), ^ll), < C, i G [0,1*], r > 0. 



Using i|4.15|) . we obtain 

|e 7 (t,r) |<C7 7 3/2 (r 



7f 



<C7,r>0,t6[0,T'] l7 6]0,l]. 



The remaining terms in H4.23[) can be estimated exactly as in the proof of Theo- 
rem 2 of ^31- But the sole change is on estimate concerning ASg/i 7 (i, x)—K(t, x) 
that involves the quantity (m 7 /r 2 — m/r 2 )(t, r). In our case, we have after dis- 
tinguishing the cases r < ro + t/^/j and r > ro + 1/\/7 : 



rn 



(t,r)<C\\ Mt)-f(t) \\ La 



C 



z 2X ^e 2 Jt,s)ds 



(4.24) 
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and using once again 1)4. 15[) . the integral that appears in (|4.24|) can be estimated 
as: 
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e 2 MM) e 2^ s)ds = g 2 7 3 J ±_ dg (J r 2 e A T ( t ,r) M ^ ( ^ ^ \ 



<C 7 3 I ' dsfe^'^Myfarfdry 
V7, 



< C 7 3 Q 7 (i)^ 6 (i) ( r + 



< c 7 . 

So, the proof of Theorem 4.2 is complete. 
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